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Abstract 

There is ample evidence of twisted magnetic structures in the solar corona. This motivates 
us to consider the magnetic twist as the cause of Alfven frequency continuum in the coronal 
loops, which can support the resonant absorption as a rapid damping mechanism for the 
observed coronal kink magnetohydrodynamic (MHD) oscillations. We model a coronal loop 
with a straight cylindrical magnetic flux tube which has constant but different densities in 
the interior and exterior regions. The magnetic field is assumed to be constant and aligned 
with the cylinder axis everywhere except a thin layer near the boundary of the flux tube 
which has an additional small magnetic field twist. Then, we investigate a number of possible 
instabilities that may arise in our model. In the thin tube thin boundary approximation, we 
derive the dispersion relation and solve it analytically to obtain the frequencies and damping 
rates of the fundamental (Z = 1) and first/second overtone (Z = 2,3) kink (m = 1) MHD 
modes. We conclude that the resonant absorption by the magnetic twist can justify the rapid 
damping of kink MHD waves observed in coronal loops. Furthermore, the magnetic twist 
in the inhomogeneous layer can cause deviations from P1/P2 = 2 and P1/P3 = 3 which are 
comparable with the observations. 
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1 Introduction 

The first identification of transverse oscillations of coronal loops was reported by Aschwanden 
et al. (1999) and Nakariakov et al. (1999) using the Transition Region and Coronal Explorer 
(TRACE) observations of 14 July 1998 in the 171 A Fe IX emission lines. Nakariakov et al. 
(1999) for a loop with length of (130 ± 6) x 10 3 km and width of (2.0 ± 0.36) x 10 3 km reported 
a spatial oscillation with period of 4.27 ± 0.13 min and decay time of 14.5 ± 2.7 min. They 
suggested resonance of the global mode as the cause of such fast damping. On 17 April 2002, 
the vertical polarization of coronal loops oscillations with period of 3.9 min and decay time of 
11.9 min were identified by Wang & Solanki (2004) using the TRACE observations in the 195 
A Fe XII emission line. According to Roberts et al. (1984), the goal of coronal seismology is to 
deduce the properties of the solar corona using observed parameters of oscillations and waves. 
For instance, Nakariakov &: Ofrnan (2001) applied a new method for determination of the local 
magnetic field strength base on the observed length, density and frequency of an oscillating 
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coronal loop. For reviews on coronal seismology, see e.g. De Moortel (2005), De Moortel & 
Nakariakov (2012), Andries et al. (2009) and Ruderman &; Erdelyi (2009). 

The theory of resonant absorption of MHD waves was first established by Ionson (1978) as 
a conceivable mechanism for heating of the solar corona. Since then, many theoretical works 
have been done to develop this theory (see e.g. Davila 1987; Sakurai et al. 1991a,b; Goossens et 
al. 1992; Steinolfson & Davila 1993). In this mechanism, energy of the global mode oscillations 
is transferred to the local Alfven perturbations within a resonance layer inside the loop. The 
necessary condition for this process is a gradient of Alfven frequency in this layer which varies 
between the interior and exterior Alfven frequencies of the loop. For a good review on the 
resonant absorption see also Goossens et al. (2011). 

Heating of the coronal loops by the resonant absorption of MHD waves was studied by Erdelyi 
& Goossens (1995). Solving the visco-resistive MHD equations of motion, they concluded that 
under coronal conditions both viscus and ohmic dissipation mechanisms are important. Erdelyi 
&; Goossens (1996) showed that the equilibrium plasma flow in coronal flux tubes affects the 
resonant absorption rate due to driving waves. 

Ruderman & Roberts (2002) investigated the resonant absorption of kink mode oscillations 
in a coronal loop with radial density inhomogeneity in the zero-beta approximation. They 
concluded that the oscillations of coronal loops are coherent only in the presence of small scale 
inhomogeneities in density. Goossens et al. (2002) showed that the damped quasi-modes give 
an accurate description of rapid damping of the observed coronal loop oscillations if the length 
scale of inhomogeneity changes from loop to loop. Also damping of quasi-modes is completely 
consistent with large estimates of the Reynolds numbers in the corona (10 14 ). 

Van Doorselaere et al. (2004) investigated the oscillations of kink modes in coronal loops 
by the LEDA numerical code (van der Holst et al. 1999). Taking into account the large 
inhomogeneity length scales, Van Doorselaere et al. (2004) showed that the rapid damping 
of oscillating coronal loops can be justified by resonant absorption without resorting to the 
Reynolds numbers smaller than the classical values. They concluded that the numerical results 
of damping rates can deviate from the analytical ones obtained in thin boundary approximation, 
by up to 25% . 

Terradas et al. (2006a) investigated the temporal evolution of resonant absorption in a one¬ 
dimensional cylindrical coronal loop. They found that when the coronal loop is excited by the 
external perturbation, the first stage of the loop oscillation has the leaky behavior. After that, 
the loop oscillates like the kink mode and then it is dissipated by the mechanism of resonant 
absorption. Terradas et al. (2006b) developed their previous work and found that considering the 
curvature of coronal loops enhances the efficiency of resonant absorption slightly. They showed 
that there are two kink modes with polarizations mostly in the horizontal/vertical directions 
with respect to the photosphere. These modes show both resonant and leaky behavior at a same 
time. 

The effect of longitudinal density stratification on the resonant absorption of MHD waves 
in coronal loops with radial density gradients has also been studied in the literature (see e.g. 
Andries et al. 2005; Karami et al. 2009). Karami et al. (2009) showed that in the zero- 
beta approximation, when the stratification parameter increases, both the period and damping 
time of the kink and fluting modes decrease but the stratification does not affect the ratio of 
frequencies to damping rates. They further showed that the ratio of fundamental period to first 
overtone one decreases from its canonical value P 1 /P 2 = 2 when the stratification parameter 
increases. 

Besides the above considerations, there are observational evidences for twisted magnetic 
fields in coronal loops. Chae et al. (2000) reported the traces of rotational motions in coronal 
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loops and suggested that the existence of an azimuthal magnetic held component that encircles 
the axis may be required to guide the rotational motions. Chae Sz Moon (2005) considered a 
model of twisted flux tube in which the force of pressure gradient is balanced by the tension force 
of the azimuthal magnetic held component. They suggested that the constriction of plasma can 
be used to determine the magnetic twist in coronal loops. Kwon Sz Chae (2008) measured the 
magnetic twist of 14 coronal loops. The magnetic twist value (f >twist of the loops was in the range 
[0.227T, 1.737r]. 

Many theoretical works have also been done on the effect of twisted magnetic helds on the 
MHD waves in coronal loops (see e.g. Bennett et al. 1999; Erdelyi Sz Carter 2006; Erdelyi Sz 
Fedun 2006, 2007, 2010; Carter Sz Erdelyi 2007, 2008; Ruderman 2007, 2015; Karami Sz Barin 
2009; Karami Sz Bahari 2010, 2012; Terradas Sz Goossens 2012). Ruderman (2007) considered 
a straight hux tube in the zero-beta approximation with a magnetic twist inside the loop pro¬ 
portional to the distance from the tube axis. Using the asymptotic analysis in the limit of 
small twists, Ruderman obtained an analytical solution for perturbations inside the loop and 
showed that the magnetic twist does not affect the standing kink modes. Karami Sz Bahari 
(2012) extended the work of Ruderman (2007) to a magnetically twisted flux tube containing 
both core and annulus regions. They showed that the frequencies and the period ratio P 1 /P 2 of 
the fundamental and first overtone nonaxisymmetric kink and fluting modes are affected by the 
twist parameter of the annulus. Terradas Sz Goossens (2012) studied the effect of magnetic twist 
on the kink oscillations of coronal loops with a piecewise parabolic twist profile. They solved 
the MHD equations using the PDE2D (Sewell 2005) code. Terradas Sz Goossens (2012) in the 
limit of small twists showed that the magnetic twist changes the polarization of the transverse 
motions of standing kink oscillation along the flux tube but does not affect its frequency. Ru¬ 
derman (2015) investigated the effect of a continues twisted magnetic held on the propagating 
kink modes in the thin tube approximation. He showed that there are two propagating kink 
waves with the same longitudinal wave numbers but opposite propagation directions, which have 
different frequencies. Ruderman (2015) called these waves “accelerated” and “decelerated” kink 
waves which have larger and smaller frequencies with respect to the well known kink frequency, 
respectively. 

Karami Sz Bahari (2010) studied the effect of twisted magnetic held on the resonant absorp¬ 
tion of MHD waves due to the radial density structuring in coronal loops. They showed that 
by increasing the twist parameter, the frequency, the damping rate and their ratio for both the 
kink and fluting modes increase. Also the magnetic twist causes the ratio of fundamental period 
to first overtone one for kink and fluting modes to be smaller than 2. 

The main goal of the present work is to study the resonant absorption of kink MHD waves 
by magnetic twist to explain the rapid damping of oscillating coronal loops and departure of the 
period ratios P1/P2 and P1/P3 from their canonical values reported by observations. To this 
aim, in section [2] we introduce the coronal loop model and find the solutions of the equations 
of motion. In section [3j we investigate a number of possible instabilities that may arise in 
our model. In section [U we use the appropriate connection formulae to obtain the dispersion 
relation. In section 0 we solve the dispersion relation, analytically. Section [6] gives the summary 
and conclusions. 
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2 Model and equations of motion 


As a simplified model for a coronal loop, we consider a straight cylindrical flux tube with length 
L and radius R. The background density profile is assumed to be 


p(r ) 


Pi, 0 < r < R, 
Pei r > R, 


( 1 ) 


where p x and p e are the constant densities of the interior and exterior regions of the flux tube, 
respectively. We define the density ratio £ = Pi/p e in the rest of the paper. 

We further assume that the background magnetic field to have a small twist in a thin layer 
and to be constant and aligned with the loop axis everywhere else, i.e. 


B = 


Bqz, 0 < r < a, 

B v (r)ip + B 0 z, a<r<R, 
B 0 z, r > R. 


Here, we set a parabolic profile for B^ as follows 


( 2 ) 


B <p(r) = Ar(r — a), (3) 

where A is a constant. Note that the jump in the value of B^ijr) across r = R implies a 
delta-function current sheet there. 

The magnetohydrostatic equilibrium equation takes the form 


d + B% 

dr ^ 2p 0 J p 0 r ’ 


(4) 


where po is the magnetic permeability of free space. Then, by integrating Eq. ([][]) and using the 
continuity of the total (magnetic plus gas) pressure at r = a and r = R, we can find the gas 
pressure as follows 


P{r) 


2m 




M_R _ ^—(r — a) 

2/i 0 P „ 12fio V ' 

x (9 r" 


— liar 2 + a 2 r + a 3 ), 


5 l 

12(io 


P- 


(3 R a — 8 aR? + 6 a?R 2 — a 


0 < r < a, 
a < r < R, 
r > R, 


(5) 


where parameter /3 = L;2 ^ is the ratio of the plasma pressure to the magnetic field pressure, 
inside the loop. 

The linearized ideal MHD equations for incompressible plasma are given by 


d5v \75P 1 

^- = -2— +-{(V x 6B) x B + (V x B) x <5B}, (6) 

dt p p 0 p 


fMB 

~df 


V x (Sv x B), 


(7) 


V • <5v = 0, 


( 8 ) 
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where <5v, <5B and 5P are the Eulerian perturbations of velocity, magnetic fields and plasma 
pressure; B and p are the background magnetic filed and the mass density, respectively. Also t-, 
4>- and ^-dependency of the perturbations are supposed to be of the form exp [i{rruj) + k z z — uit)] 
where k z = ^ is the longitudinal wave number. Here m and l are the azimuthal and longitudinal 
mode numbers, respectively, and w is the mode frequency. So the perturbations are of the form 


5v(r, 4>, z, t ) = Sv(r) exp [i(m(f> + k z z — cut)], 
<5B(r, i f>, z, t) = <5B(r) exp [i(m<j) + k z z — ut)], 
(5P(r, (j>, z, t ) = <5P(r) exp[i(m</> + k z z — ut)\. 


(9) 


Following Bennett et al. (1999), substituting the perturbations @ into Eqs. dU|) and Q and 
doing some algebra we get 


dr 2 P 1 rD dr 


rD 


dr 


5p 


+ < (§) + w {C2 ° 3 ~ c2) 

p D d Ci 
t r = C- 3 Tr 6p+ C- 3 6p ’ 


(10a) 

(10b) 


where dp = 5P + B • 6 B//xo and = —5v r /iu are the Eulerian perturbation of total pressure 
and the Lagrangian displacement in the radial direction, respectively and 


D = p(uj 2 -u >|), 
Ci = - 


2 mB m pm „ , „ 

- 2~ (—+ k z B 2 

por z \ r 


c 2 = - 


m 


+ k 


( 11 ) 


C 3 = D 2 + D^b- 

po dr 

Here, the Alfven frequency, u A , is defined as 

1 


B 


4 B 2 

2 

-2 P^A- 

Por 2 


u A (r) = -=(-B v + k z B z ). 
yfpPp \r ) 

Putting Eqs. m-m into m gives the profile of Alfven frequency as follows 


( 12 ) 


u A (r) = < 


Bnk z 
VMOPi 1 


0 < r < a, 


= (mA(r - a) + k z B^j , a < r < R, 

r > R. 


V^OPi 
Bp k, 

VMOPe ’ 


(13) 


Here we use the twist parameter dehned as a = B ' p ^ ■ For a special value of a = a c , the Alfven 
frequency is continues at the tube boundary (r = R). If a ^ a c there would be a gap in the 
Alfven frequency profile across the boundary. It is straightforward from Eq. (1131) to show that 


Rk, 


Qfr — 


m 


(A 


- 1 


(14) 
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Figure 1: Azimuthal component of the background magnetic held, Alfven frequency of the 
fundamental (1 = 1) kink (m = 1) mode, background Alfven speed and background plasma 
pressure versus the fractional radius r/R. The loop parameters are: L = 10 5 km, R/L = 0.01, 
a/R = 0.95, ( = 2, p\ = 2 x 10 -14 g cm -3 , j3 = 0.1 and Bq = 100 G. Here a = a c = 0.013 and 
va { = Bo/^/^opi = 2000 km s -1 . 


Figured] shows the azimuthal component of the background magnetic held, the Alfven frequency 
of the fundamental (l = 1) kink (m = 1) mode, the background Alfven speed, va(t) = B/^/pop, 
and the background plasma pressure for the twist parameter a = a c = 0.013. Note that although 
for a > a c , there is a sudden drop in the Alfven frequency ua( r) at the tube boundary, the 
interior Alfven speed va(t) in the present model is still smaller than the exterior one which is 
likely to occur under coronal conditions. 

From Eqs. (H0aD - (ll2D . it is clear that when oj 2 = i.e. D = 0, the equations of motion 
in the inhomogeneous layer (a < r < R) become singular in the presence of magnetic twist. 
Therefore, resonant absorption can occur not only by the radial density inhomogeneity (like 
previous works, see e.g. Ruderman Sz Roberts 2002; Kararni et al. 2009; Karami Sz Bahari 2010; 
Ruderman 2011; Ruderman Sz Terradas 2013; Soler Sz Terradas 2015) but also by the magnetic 
twist (present work). 

In the untwisted regions, 0 < r < a and r > R, Eqs. (llOal) and (llObD are reduced to the 
following equations 


d 2 Sp 
dr 2 


1 d 5p 

~ l - 

r dr 


m 


k z + — ) 5p = 0 , 


(15a) 


1 d 5p 
p(u 2 — lo\) dr ' 


(15b) 


Solutions for equations (ll5aD and (115bl) in the interior (0 < r < a) and exterior (r > R) regions 
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are obtained as 


6p(r) = 


Ail m (k z r), 0 < r < a, 
A e K m (k z r), r > R, 


(16a) 




A e Pe( J-^ Ae ) K Ukz r ), T > R. 


0 < r < a 


(16b) 


Here I m and K m are the modified Bessel functions of the first and second kind, respectively. 
Also “ / ” on I m and K m represents a derivative with respect to their arguments. The constants 
Ai and A e are determined by the appropriate boundary conditions. 

3 Stability constraints on the model 

In what follows, we are interested in investigating a number of possible instabilities that may 
arise in our model. 

3.1 Kink instability 

If the magnetic twist value in the loop, which is defined as </>twist = = 27rIVtwist exceeds 

a critical value < f> c , then the loop becomes kink unstable (see e.g. Shafranov 1957; Kruskal et 
al. 1958). Here, A r tw i s t is the number of twist turns in the loop. For force-free magnetic fields 
of uniform twist, a critical twist of <^>twist > 3.3-7T (or 1.65 turns) was found to lead to kink 
instability (Hood & Priest 1979), while the critical value ranges between 27 t and 6ir for other 
types of magnetic fields (see e.g. Aschwanden 2005; Priest 2014). Numerical MHD simulations 
of an increasingly twisted loop system demonstrated linear instability of the ideal MHD kink 
mode for twist angles in excess of ~ 4.87T (or 2.4 turns) (Mikic et al. 1990). Baty & Heyvaerts 
(1996) examined the kink instability of a radially localized twist profile and obtained <f> c = 57t. 
In our model to avoid the kink instability, following Hood & Priest (1979) we consider cj) c = 3.37T 
(or 1.65 turns). This yields an upper limit for the twist parameter as a max = ( R/L)(j) c . 

3.2 Kelvin-Helmholtz instability 

The Kelvin-Helmholtz instability (KHI) can occur during the kink oscillations of coronal loops. 
The kink wave, whose initial energy is mostly transverse, converts into azimuthal waves locally 
resembling torsional Alfven waves, which are finely localized around the tubes boundary layer 
(Verth et al. 2010; Arregui et al. 2011). Azimuthal motions are thus amplified and introduce 
velocity shear, making them prone to be unstable to the KHI (Heyvaerts & Priest 1983; Soler 
et al. 2010). Terradas et al. (2008) using numerical simulations showed that azimuthal shear 
motions generated at the loop boundary during kink oscillations can give rise to a KHI, but 
this phenomenon has not been observed to date. The KHI can also extract the energy from the 
resonant layer and convert it into heat through viscous and ohmic dissipations at the generated 
vortices and current sheets (Antolin et al. 2015). Soler et al. (2010) pointed out that the 
presence of a small azimuthal component of the magnetic field can suppress the KHI in a stable 
coronal loop. They showed that the required twist is small enough to prevent the development 
of the pinch instability. A weak twist of magnetic field lines is very likely and realistic in the 
context of coronal loops. 
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3.3 Resistive kink instability 

The effect of resistive diffusion on the ideal kink mode yields a so-called resistive kink instability 
which is a reconnecting process. The resistive kink becomes important when the loop is twisted 
too much (see e.g. Biskamp 2000; Wesson 2004; Priest 2014). From Eq. ([7]). it is possible to 
have 

Vx(ivxB) = (B-V)5v - B(V • <Sv) = 0, (17) 

at a specific location r = r s . In this case, the diffusion term t]\7 2 5'B neglected from the right hand 
side of Eq. (0 becomes important and consequently the field lines diffuse through the plasma 
and reconnect. Here rj is the magnetic diffusivity. For an incompressible plasma, inserting 
Eq. 0 into (1171) gives rise to k • B| = 0, where k and B are the wave vector and the 
background magnetic field, respectively. Using Eq. 0, the necessary condition for the resistive 
kink instability takes the form 

k • B| = + k z B z (r s ) = 0. (18) 

l / —/ s rp 

1 s 

In Fig. El we plot k • B for the fundamental (/ = 1) kink (m = 1) mode. Figure illustrates that 
the term k • B cannot be zero and consequently the resistive kink instability cannot occur in our 
model. 


3.4 Tearing mode instability 

The tearing mode instability can occur in a thin current sheet when the driving force of the 
inflow exceeds the opposing Lorentz force (see e.g. Furth et al. 1963; Goldstone &; Rutherford 
1995; Magara & Shibata 1999; Tenerani et al. 2015). As a consequence of non-zero resistivity, 
the magnetic field lines tear and reconnect in the current sheet. According to Furth et al. (1973), 
the smallest growth time of the tearing mode in the cylindrical flux tube is given by about 


2/5 3/5 
T t~ T A T d > 


(19) 


where ta = I s /va ; and = l 2 /r] are the Alfven and magnetic diffusion time scales, respectively. 
Here l s is a half-thickness of the current sheet and va { = -Bo/\/PoPi is the interior Alfven speed. 
In our model, due to having a rotational discontinuity of the background magnetic field at the 
tube boundary, we have a delta-function current sheet. Therefore, from Eq. m when Zg —x 0 
the growth time of the tearing process goes to zero. It should be noted that the delta-function 
current sheet in our model is used to approximate a finite width current layer. In reality, the 
thickness of such current sheets in the solar corona takes place in the range of macroscopic 
values (10 km, for example), see e.g. Magara & Shibata (1999). For the present model, if 
we take T = 5 x 10 6 K, q = 10 9 T -3 / 2 = 0.09 m 2 s” 1 (Priest 2014), = 2000 km s" 1 

then for 21q ~ (1 — 10) km (Magara &; Shibata 1999) we estimate the tearing growth time as 
Tt/Pkmk ~ (3 — 123). Here -Pkink = 27r/u;kmk — 87 s is the period of the fundamental kink mode 
oscillation (see Eq. 1361) . Therefore, by choosing an appropriate thickness for the current sheet, 
the tearing mode instability can be avoided in our model during the resonant absorption of the 
kink oscillations. 

It is worth noting that considering a continuous background magnetic field like the profile 
used in the work of Hood et al. (2016), allows us to have a twist in the tube without having 
the delta-function current sheet. This will make it easier to avoid the growth of the tearing 
instability. However, for a certain density distribution (like the one we use in this paper) the 
model of Hood et al. (2016) does not yield a monotonic function of the background Alfven 
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Figure 2: Variation of k • B for the fundamental (l = 1) kink (m = 1) mode versus the fractional 
radius r/R. Here a = 0.1 and the other auxiliary parameters as in Fig. [0 


frequency. As a result it may give rise to two resonant layers. In this case, in order to use the 
connection formulae (see section 4) at the resonance layers, we need to solve the MHD equations 
in the twisted regions around the resonance points. However, this is beyond the scope of this 
paper and we leave it for future work. 


4 Connection formulae and dispersion relation 


Here, we do not solve Eq. (llOal) in the inhomogeneous layer (a < r < R), where the singularity 
occurs due to the existence of the magnetic twist. Instead, the solutions inside and outside of 
the tube can be related to each other by the connection formulae introduced by Sakurai et al. 
(1991a). To check the validity of the connection formulae, the radius Sa of the region around 
the resonance point that connects the solutions of the perturbations of the interior and exterior 
of the flux tube, must obey the following condition (see Goossens et al. 2011) 


S A <^h = 


& 


r=r A 


( 20 ) 


In our work, we set S A = In Fig. [3l we plot S A and h versus the twist parameter 

a = — Bo' Fig ure clears that for the range of magnetic twist considered in this study, the 
condition (12011 holds. Therefore, in the thin boundary approximation we can use the connection 
formulae around the inhomogeneous layer. According to Sakurai et al. (1991a) the jumps across 
the boundary (resonance layer) for 5p and are 


where 


[5p] = 

Kr] = 


in 2 B ip (r A )Bz{r A )f B {rA) 


m 


Mo p\r A B 2 (r A ) 
3B{TA) CUM, 


C A (r A ), 


|A| p[B 2 (r A ) 


ri x f \ VbB v B z / x 

C A = 9B0p(r) --— £r(r), 

P 0 r A 
7 Tl 

Ib — Bp + k z B z , 
r 

Jfl 

9 b — B z k z I3 ( p , 

r 

A = _ ^ (r) ’ 


(21a) 

(21b) 


( 22 ) 


9 













10" 1 
10- 2 
10" 3 
10- 4 
10" 5 

0.02 0.04 0.06 0.08 0.10 

a 

Figure 3: Variations of h. S A = 0.025 and thickness of the resonance layer 5 A versus the twist 
parameter a for the fundamental (/ = 1) kink (m = 1) modes. Here h, 5 a and S A are in units 
of the loop radius R = 1000 km. Auxiliary parameters as in Fig. [0 



and rA is the location of the resonance point. Notice that the ideal MHD solutions for the 
untwisted regions, i.e. inside (0 < r < a) and outside (r > R) the flux tube, are valid in the 
jump conditions (I21af) and (121bl) only if the following condition is established (see Stenuit et al. 
1998; Goossens et al. 2011) 

(23) 

where 5a = |jSj( z/ + v) )| 1//3 is a measure of the thickness of the 

resonance layer (see Sakurai et al. 1991a). Here v = p = ^— and 1Z V are the kinematic 
viscosity, magnetic diffusivity and viscus Reynolds number, respectively. The classical values of 
viscus and resistive Reynolds numbers of the solar corona are about 10 14 and 10 13 , respectively 
(see Colub & Pasachoff 1997). In Fig. [3l we also plot 5a for the aforementioned values of the 
Reynolds and Lundquist numbers. Figure clears that the condition (fTTl) is respected. 

Substituting the ideal solutions (llfia.l) and (I16bli in the jump conditions (121af) and (121bl) and 
eliminating A[ and A e . one can find the dispersion relation as 

do(oj) + di(u) = 0, (24) 


where 


do{u>) 


k z K' m (k z R ) 

p e (uj 2 - K m (k z R ) 

kz I m (k z a) 

Pi(u 2 -uj 2 a .) I m {k z a)' 


(25) 


and 


di(u>) = 
9b 


1 7T 

]A| 


k z 


2 f B B v K' m (k z R ) 


9B 


p x Bl p e (uj 2 - p, 0 pir A B 0 K m (k~R) 
I m {k z r A ) k z 2f B B v B 0 1' m (k z r A ) 


i(k z a) pi(a ; 2 - w\.) p 0 r A I m (k z a ) 


(26) 


Here u) = oj — i'y, in which a; and 7 are the mode frequency and the corresponding damping rate, 
respectively. 
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Using Eq. (fT3l) in the last relation of Eq. (1221) we get 


Bq 2 ma 

MoPi R(R ~ a) 


^ma 


r — a 
R(R — a) 



(27) 


which shows that in the limit of a —>• R, i.e. when the twisted annulus region is removed, then 
A goes to infinity. As a result, in Eq. (1261) . d\ = 0. In this case there is no resonant absorption 
and the dispersion relation (1241) is the same as Eq. (25) in Bennett et al. (1999). 

In the limit of thin tube approximation, using the first order asymptotic expansions for 7 m (x) 
and K m (y) one can get 


Im(k z r A ) _ / 7A \ m 
ImM ~ \ a J 

I' m {k z r A ) _ m / 7 a \ m 

Im{k z a) ~ k z r A U/ 

dm (k z ci) k z a 
K'm{k z R) ^ m 
K m (k z R) k z R 


(28) 


Using the above approximations, the dispersion relation (1241) reduces to 


Pi (oj 2 - uj 2 a .) p e (uj 2 - uj 2 a J 
R a 



x 




2 N , 2mf B B lf \ 

+ I^mr) 

2f B B^B 0 \ _ 

kmr'i ) 


(29) 


Notice that we look for the frequencies in the range of uj Ai < uj < uj Ae , in which the resonant 
absorption occurs. Here uj a . = k z &= and uj a = k z —M^= are the interior and exterior Alfven 

1 z y/fiOPi Sie z y/POpe _ 

frequencies, respectively. In the next section, we solve Eq. (1291) analytically to obtain the 
frequency and damping rate of the kink MHD modes. 


5 Analytical results 


Here, we are interested in studying the effect of the twist parameter a on the frequencies, uj, and 
the damping rates, 7 , of the kink (m = 1) MHD modes in a resonantly damped coronal loop. 
To this aim, we need to solve Eq. (1291) . First, we use the dimensionless quantities f = r/R, 
L = L/R, B = B/Bq and uj = uj/(v a JR). Thus, Eq. (1291) can be recast in dimensionless form 
as 


((uj 2 


2 \ 1 1 r~2 2 \ 

~ w A i ) + “( a; ~ bJ AJ 


in ( r A 


I Al 


a 


( 9B(rf - U 2 Ae ) + 


2(mf B B y 
r A 




2 /bM =q 
r A J 


(30) 
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where we have dropped the bars for simplicity. Here q = a/R. Equation (1301) yields a quadratic 
equation for Co 2 as follows 


ciio 4 + c 2 £u 2 + c 3 = 0 , 


(31) 


where 


ci = - 


, . m 2 

m r A \ g% 


|A| \ q ) m ’ 


TT ( W ii + W Ae) + 


X 


2 


m 


(l - C^a) 


l 


ITT / TA 


c 3 - -c^ ( 1 + - J - ]K[ VV ) £ 


9% 4 . 2 fB9BB<pU A . 


X 


- u A- 

m 1 


+ 


-(1 - r A ) 


'A 'A 

It is straightforward to find the solutions of Eq. m, Cj = u — * 7 , as follows 

1 / /---— \V 2 


w± = 


9l| + 3 2 ± 


7± = 


(\/x± 2 + 0± + x±) 




where 


X=t C 3 R 91 -l- H 03 ^ 3 -, 
©± = C3 : 91-t - c 3r 3±, 


= -c 2r ± 


\J c 4 + c 5 + C 4 


1/2 


3± = — c 2 


± ^ \Z C 4 + C 5 ~C 4 ^j 7 


with 


c 4 — c 2 r — c 2i + 4ci r C 3l , 
C 5 = 2c 2r c 2i — 4 cijC3 R . 


(32) 


(33) 


(34) 


(35) 


Here, the subscripts R and / denote the real and imaginary parts of a given complex quantity. 

For a typical coronal loop, we take L = 10 5 km, R/L = 0.005,0.01, C = Pi/Pe = 2,4, 
Pi = 2 x 10 ~ 14 g cm ~ 3 and Bo = 100 G. So, one finds v Ai = Bo/^/popi = 2000 km s -1 in this 
loop. We assume that the magnetic twist takes place in a thin layer of thickness R — a = 0.05 R. 
So we have q = 0.95. In the thin boundary approximation, we can set the location of the 
resonance point at the tube surface, i.e. r A ~ R. 
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Using Eq. (1331) we obtain two roots of oj± which their values take place in the range of 
ojAi < W- < uja b and cj + > loa b - Therefore, cj_ is our physical root and should be ruled 
out, because it does not give rise to the resonant absorption. As we have already mentioned, 
in the limit of q —>• 1 (i.e. a —>• R ), we have A —>• oo. In this case, from Eq. (1321) we obtain 
ci = C 2 : = C 3 j = 0. Hence, in the absence of twisted annulus region, form Eqs. P51) - (1551) we get 


0J- = 



k-W — k4dnk> 


7 - = 0 , 


(36) 


where Wkink is the kink mode frequency of an untwisted loop in the thin tube approximation 
(see e.g. Goossens et al. 2009). 

Here, we are interested to obtain a minimum value for the twist parameter a = a m ; n which 
is required to excite the resonance. The exact value of a m i n is determined as the root of the 
following equation 


W-\ a=ainin — U>A{R)\a=a min - (37) 

We solve this numerically and conclude that the result obtained for a m ; n satisfies the following 
relation 


|ct=ct m i n — ^kink- (38) 

Equating right hand sides of Eqs. (EZD and (1381) , one can get an approximate expression for a m i n 
as 



For instance, taking £ = 2 and R/L = 0.01, for the fundamental (/ = 1), first overtone (l = 2) 
and second overtone (l = 3) kink (m = 1) modes, Eq. (139)) gives a m i n =0.0049, 0.0097 and 
0.0146, respectively. The relative error of the results of Eq. (1391) with respect to the exact 
solutions of Eq. (1371) is about 0.6%. 

Figures HI [5l and [ 6 ] illustrate the frequencies (w_ = co m i), the damping rates (y_ = 7 m i) 
and the ratio of the frequency to the damping rate ( 0Jmi/lmi ) of the fundamental (1 = 1) and 
first/second overtone (l = 2,3) kink (m = 1) modes versus the twist parameter a, respectively. 
Figures show that (i) the frequencies and damping rates increase when the twist parameter 
increases, (ii) The ratio of the oscillation frequency to the damping rate w/7 decreases when the 
twist parameter increases. The result, interestingly enough, is that for the fundamental (l = 1) 
kind (m = 1) mode, for the twist parameter a = 0.0147 we obtain wn/(27r7n) = 3 which is in 
good agreement with the observations reported by Nakariakov et al. (1999), Wang & Solanki 
(2004), and Verwichte et al. (2004). Here, the behaviour of u and 7 versus the twist parameter 
a are the same as that obtained by Karami & Bahari (2010) but for the resonant absorption 
due to the radial density structuring, (iii) For a given a and (, when R/L increases then co 
increases, 7 decreases and u/ 7 increases, (iv) For a given a and R/L, when ( increases, the 
frequency and the ratio u/ 7 increase. However, this behavior for the damping rate holds only 
for large values of the twist parameter. Moreover, increasing ( does not affect the ratio ui /7 
for large values of a. (v) Note that based on Eq. (1391) . a m ; n is a function of m, l, £ and R/L. 
Therefore, for each set of these parameters, the start points of the diagrams in Figs. [7] to [U] have 
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Figure 4: Frequency of the fundamental (l = 1) kink (m = 1) mode (top), its damping rate 
(middle) and the ratio of the oscillation frequency to the damping rate (Bottom) versus the twist 
parameter a. Solid line: £ = 2 ,R/L = 0.01; Dashed line: ( = A,R/L = 0.01; Dash-Dotted line: 
( = 2 ,R/L = 0.005. Other auxiliary parameters as in Fig. [TJ Both frequencies and damping 
rates are in units of vaJR = 2 rad s —1 . 


different a m ; n values. Also to avoid the kink instability in our model, following Hood & Priest 
(1979) we consider the twist value 0twist < 4> c = 3 . 37 T and consequently obtain an upper limit 
for the twist parameter as a max = ( R/L)(f> c . Hence, the diagrams in Figs. [4] to [HI with different 
R/L have different a max cut-offs. 

Figure [7] presents variations of uj, 7 and uj/ 7 for the fundamental (l = 1) kink (m = 1) mode 
versus the thickness of the twisted layer, d = R — a. Figure shows that when d increases both 
frequency and damping rate increase but the ratio of the frequency to the damping rate decreases. 
This behaviour holds also for the first and second overtone (l = 2,3) kink modes. Figure [7] 

clarifies that in the limit of d —> 0 (i.e. a —>• R), we have wn —> Wkmk = \J— 0.03628 

(in units of 2 rad s -1 ) and 711 —>• 0. Note that the result of cj /7 for the fundamental mode of 
kink oscillations is in good agreement with that obtained for the resonant absorbtion due to the 
radial density structuring (see e.g. Ruderman & Roberts 2002; Ruderman & Terradas 2013). 

In Fig. [H the period ratios of the fundamental to the first overtone, P1/P2, and to the second 
overtone, P1/P3, of the kink (m = 1) modes are plotted versus the twist parameter a. Figure 
shows that (i) when the twist parameter increases, the values of P1/P2 and P1/P3 decrease from 
their canonical values, 2 and 3, respectively, (ii) The density ratio ( = Pi/p e does not affect the 
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Figure 5: Same as Fig. 3] but for the first overtone (Z = 2) kink (m = 1) mode. 
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Figure 6: Same as Fig. [5] but for the second overtone (/ = 3) kink (rri = 1) mode. 
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Figure 7: Frequency of the fundamental (l = 1) kink (m = 1) mode (top), its damping rate 
(middle) and the ratio of the oscillation frequency to the damping rate (bottom) for a = 0.05 
versus the thickness of the twisted layer, d (in units of R = 1000 km). Other auxiliary parameters 
as in Fig. |T] 
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Figure 8 : Period ratios of the fundamental to the first overtone, P 1 /P 2 , and to the second 
overtone, P 1 /P 3 , of the kink (m = 1 ) modes versus the twist parameter a. Solid line: ( = 
2,R/L = 0.01. Dashed line (which overlaps the solid one): £ = 4 ,R/L = 0.01. Dash-Dotted 
line: ( = 2 ,R/L = 0.005. Other auxiliary parameters as in Fig. [TJ 




d 


Figure 9: Period ratios of the fundamental to the first overtone, P 1 /P 2 , and to the second 
overtone, P\ /P 3 , of the kink (m = 1) modes for a = 0.05 versus the thickness of the twisted 
layer, d, (in units of R = 1000 km). Other auxiliary parameters as in Fig. [TJ 
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Table 1: Magnetic twist parameter a = B^/Bz and corresponding twist value </>twist = {ji) a 
predicted by our model for some observational period ratios P1/P2 and P1/P3 of the kink (m = 1 ) 
modes. Other auxiliary parameters as in Fig. [TJ 


Reference 

P 1 /P 2 

a 

0twist / 7T 

Van Doorsselaere et al. (2007) 

1.795 ±0.051 

0.0718 

2.28 

Van Doorsselaere et al. (2009) 

1.980 ±0.002 

0.0295 

0.94 

Ballai et al. (2011) 

1.82 ±0.02 

0.0678 

2.16 

Reference 

P 1 /P 3 

a 


Van Doorsselaere et al. (2009) 

2.89 ±0.14 

0.0460 

1.46 


period ratio. Notice that the solid curve overlaps the dashed one. (iii) For a given a and (, the 
period ratio decreases when R/L decreases. In Fig. [21 we plot variations of the period ratios 
P1/P2 and P1/P3 versus the thickness of the twisted layer d. Figure [9] shows that the period 
ratios decrease with increasing d. 

Note that the deviations of the period ratios P1/P2 and P1/P3 from their canonical values 
provide useful information about the magnetic twist structuring within the loop. For instance, 
the observed value P1/P2 = 1.795 reported by Van Doorsselaere et al. (2007) can be justified with 
the twist parameter a = 0.0718. Van Doorsselaere et al. (2007) also reported the identification 
of the second overtone of the kink mode by re-analyzing the transverse oscillations of coronal 
loops observed by TRACE on 13 May 2001. The period ratio of fundamental to second overtone 
kink mode found to be P1/P3 = 2.89 which can be justified with the twist parameter a = 0.0460. 
Table Q] summarizes the twist parameters and the corresponding twist values, ^twist = )a, 

predicted by our model in order to justify some observed period ratios P1/P2 and P1/P3 of the 
kink (m = 1) modes. Note that the magnetic twist values 4>t w ist predicted by our model are in 
the range of the observational values reported by Kwon & Chae (2008) and Wang et al. (2015). 

6 Conclusions 

Here, we investigated the resonant absorption of kink MHD modes by magnetic twist in coronal 
loops. To this aim, we considered a thin straight cylindrical flux tube with a twisted magnetic 
field in a thin layer at the boundary of the loop and a straight magnetic field everywhere else. 
The magnetic twist causes a radial Alfven frequency gradient and consequently gives rise to 
the resonant absorption. We assumed the plasma density to be constant but different in the 
interior and exterior regions of the loop. We obtained the solutions of ideal MHD equations 
for the interior and exterior regions of the tube. Then, we investigated a number of possible 
instabilities that may arise in our model and concluded that these instabilities can be avoided in 
the present model. We also derived the dispersion relation by using the appropriate connection 
formula introduced by Sakurai et al. (1991a). In thin tube thin boundary approximation, we 
solved analytically the dispersion relation and obtained the frequencies and damping rates of 
the fundamental (l = 1) and first/second overtone (l = 2,3) kink (m = 1) modes. Our results 
show the following. 

• The frequencies and damping rates increase when the twist parameter increases. 

• The ratio of the fundamental frequency to its corresponding damping rate of the kink 
(m = 1) modes can well justify the rapid damping of kink MHD waves (w/( 27 r 7 ) ~ 3) 
reported by the observations. This confirms the high efficiency of resonant absorbtion due 
to the magnetic twist. 
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• For a given twist parameter a and density ratio £ = p\/p e by increasing R/L, the frequen¬ 
cies increase, the damping rates decrease and the ratio of the frequency to the damping 
rate increases. 

• By increasing the thickness of the twisted layer, the frequency and the damping rate 
increase but the ratio of the frequency to the damping rate and the period ratio decrease. 

• For a given a and £, the period ratio decreases with decreasing R/L. Furthermore, the 
density ratio £ does not affect the period ratio. 

• The period ratios P 1 /P 2 and P 1 /P 3 with increasing the twist parameter, decrease from 2 
and 3, respectively. For some special values of the twist parameter, the values of P 1 /P 2 
and P 1 /P 3 predicted by our model can justify the observations. 
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